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ABSTRACT 

In this paper we examine an interesting connection between the generalized Volterra lattices of 
Bogoyavlensky and a special case of an integrable system defined by Sklyanin. The Sklyanin 
system happens to be one of the cases in the classification of Kozlov and Treshchev of Birkhoff 
integrable Hamiltonian systems. Using this connection we demonstrate the integrability of 
the system and define a new Lax pair representation. In addition, we comment on the bi- 
Hamiltonian structure of the system. 

PACS numbers 02.301k, 02.30Hq. 

Keywords: Bogoyavlensky- Volterra lattices, Birkhoff integrable systems, Toda lattice, 
Poisson brackets. 



1 Introduction— Birkhoff Integrable systems 



In this paper we examine an interesting connection between the generalized Volterra lattices 
of Bogoyavlensky [6] and a special case of an integrable system defined by Sklyanin [27]. 
The Sklyanin system happens to be one of the systems in the classification of Kozlov and 
Treshchev of Birkhoff integrable Hamiltonian systems. Using this connection we are able to 
prove the integrability of the system and define a new Lax pair representation different from 
the one in [26]. In addition, we comment on the multi-Hamiltonian structure of the system. 
This connection was discovered in an effort to connect the Volterra D n system with the 
corresponding Toda D n system as in the case of the other classical Lie groups. In contrast to 
the other cases, the Volterra D n system does not correspond to the Toda D n system under 
the procedure of Moser but it actually corresponds to a special case of the Sklyanin system 
mentioned above. A Miura type transformation between the Volterra and Toda D n systems 
is therefore still an open problem. 

We begin with the following more general definition which involves systems with exponential 
interaction: Consider a Hamiltonian of the form 

ff4(p,p) + E« (v " q) - (i) 

z i=l 

where q = (qi, ■ ■ ■ ,q n ), P — (pi, • • ■ ,Pn), v i, • • • , vn are vectors in R n and ( , ) is the 
standard inner product in R n . The set of vectors A = {vi, . . . ,v N } is called the spectrum 
of the system. 

Let M be the N x N matrix whose elements are 

Mij = (vi,Vj) . 

Hamilton's equations of motion can be transformed by a generalized Flaschka transformation 
to a polynomial system of 2N differential equations. The transformation is defined as follows: 

Oi = -e^, 6i = (v i>P ). (2) 
We end-up with a system of polynomial differential equations 

«fc = akbk 

N 

b k = M ki a i ■ (3) 
i=i 

Equations (3) admit the following two integrals 

N N 

F 1 = ]TAA, F 2 = n^ (4) 

i=i i=i 

provided that there exist constants A« such that X^I=i K^i = 0. Such integrals always exist 
for N > n. One can define a canonical bracket on the space of variables (a,i,bi) by the 
formula 

{bi,aj} = (vi.VjOaj 

and all other brackets equal to zero. The integrals F\ and Fi are Casimirs of system (3). 
An interesting case of (1) occurs when the spectrum is a system of simple roots for a simple 
Lie algebra Q. In this case N = I = rankC?. It is worth mentioning that the case where 



N, n are arbitrary is an open and unexplored area of research. The main exception is the 
work of Kozlov and Treshchev [19] where a classification of system (1) is performed under 
the assumption that the system possesses n polynomial (in the momenta) integrals. We 
also note the papers by Ranada [26], Annamalai, Tamizhmani [3] and Emelyanov [10]. Such 
systems are called Birkhoff integrable. For each Hamiltonian in (1) we associate a Dynkin 
type diagram as follows: It is a graph whose vertices correspond to the elements of A. Each 
pair of vertices Vj, Vj are connected by 

4(v„v J ) 2 

(Vi,Vi)(Vj,Vj) 

edges. 

Example: The classical Toda lattice corresponds to a Lie algebra of type A n ^\. In other 
words iV = / = n — 1 and we choose A to be the set: 

Vl = (1,-1,0,...,0), ,v n _i = (0,0,..., 0,1,-1) . 

The graph is the usual Dynkin diagram of a Lie algebra of type A n _i. The Hamiltonian 
becomes: 

H( qi , ...,q n , Pl , ... , Pn ) = £ \p\ + ]T e*"** 1 , (5) 

i=l 1 i=l 

which is the well-known classical, non-periodic Toda lattice. This system was investigated 
in [12], [13], [15], [21], [23], [28] and numerous other papers that are impossible to list here. 
This type of Hamiltonian was discovered by Morikazu Toda [28] . The original Toda lattice 
can be viewed as a discrete version of the Korteweg-de Vries equation. It is called a lattice 
as in atomic lattice since interatomic interaction was studied. This system appears also in 
Cosmology, in the work of Seiberg and Witten on supersymmetric Yang-Mills theories and 
it has applications in analog computing and numerical computation of eigenvalues. But the 
Toda lattice is mainly a theoretical mathematical model which is important due to the rich 
mathematical structure encoded in it. The Toda lattice is integrable in the sense of Liouville. 
There exist n independent integrals of motion in involution. These integrals are polynomial 
in the momenta. 

As we mentioned earlier, the Toda lattice was generalized to the case where the spectrum 
corresponds to a root space of an arbitrary simple Lie group. These systems generalize the 
usual finite, non-periodic Toda lattice (which corresponds to a root system of type A n ). This 
generalization is due to Bogoyavlensky [4]. These systems were studied extensively in [17] 
where the solution of the systems was connected intimately with the representation theory 
of simple Lie groups. There are also studies by Olshanetsky and Perelomov [24] and Adler, 
van Moerbeke [1]. 

It is more convenient to work, instead with the space of the natural (q, p) variables, with the 
Flaschka variables (a, b) which are defined by: 

a, = |e^) i = l,2,...,N (g) 
bi = -\pi i = 1,2, . . .,n . 

We end-up with a new set of polynomial equations in the variables (a,b). One can write 
the equations in Lax pair form, see for example [25]. The Lax pair (L(t), B(t)) in Q can be 
described in terms of the root system as follows: 



I I 

L 0) = S kit)* 1 * + E Oi(*)( e Oi + e -aJ , 

1=1 1=1 
i 

B (t) = £a»(*)( e ai _ e -"J • 
i=l 

As usual /i aj is an element of a fixed Cartan subalgebra and e ai is a root vector corresponding 
to the simple root ctj. The Chevalley invariants of Q provide for the constants of motion. 
The first important result in the search for integrable cases of system (1) is due to Adler 
and van Moerbeke [2]. They considered the special case where the number of elements in 
the spectrum A is n + 1 (i.e., N — n + 1). Furthermore, they made the assumption that any 
n vectors in the spectrum are independent. Under these conditions a criterion for algebraic 
integrability is that 

2 ( v i. v j) m 
(Vi, v,) 1 ) 

should be in the set {0, — 1, —2, . . .}. The classification obtained corresponds to the simple 
roots of graded Kac-Moody algebras. The associated systems are the periodic Toda lattices 
of Bogoyavlensky [4]. The complete integrability of these systems using Lax pairs with a 
spectral parameter was already established in [1]. The method of proof in [2] is based on 
the classical method of Kovalevskaya. 

Sklyanin [27] pointed out another integrable generalization of the Toda lattice: 

n n— 1 

H{q u . . . , q n , Pl , . . . , Pn ) = o P 2 i + E e9l ~ q ' +1 + «i e91 + + a ne~ qn + /?„e" 2 ^ . (8) 

i=i 1 i=i 

He obtained this system by means of the quantum inverse scattering R-matrix method. This 
potential will be the main focus of this paper. 

The next development in the study of system (1) is the work of Kozlov and Treshchev on 
Birkhoff integrable systems. A system of the form (1) is called Birkhoff integrable if it has 
n integrals, polynomial in the momenta with coefficients of the form 

whose leading homogeneous forms are almost everywhere independent. We remark that in 
the definition given in the book of Kozlov [20] there is no assumption on involutivity of the 
integrals. In [19] it is proved that the polynomial integrals are in involution. The terminology 
has its origin in the work of Birkhoff who studied the conditions for the existence of linear 
and quadratic integrals of general Hamiltonians in two degrees of freedom. A vector in A is 
called maximal if it has the greatest possible length among all the vectors in the spectrum 
having the same direction. Kozlov and Treshchev proved the following theorem: 

Theorem 1 Assume that the Hamiltonian (1) is Birkhoff integrable. Let Vj be a maximal 
vector in A and assume that the vector Vj G A is linearly independent of v. Then 

2(vi, Vj) 

(Vi, 

lies in the set {0, —1, —2, . . .}. 



Note that the condition of the theorem is exactly the same as condition (7) of Adler and van 
Moerbeke. Of course theorem 1 is more general since there is no restriction on the integer N 
(the number of summands in the potential of (1)). It turns out, however, that N cannot be 
much bigger than n. In fact, it follows from the classification that N < n + 3. A system of 
the form (1) is called complete if there exist no vector v such that the set Au{v} satisfies the 
assumptions of theorem 1. In [19] there is a complete classification of all possible Birkhoff 
integrable systems based on theorem 1. The Dynkin type diagram of a complete, irreducible, 
Birkhoff integrable Hamiltonian system is isomorphic to one of the following diagrams: 




0) 



Remark 1 In the list of diagrams we have omitted some cases that occur as sub-graphs of 
diagrams (a)-(k) (by truncating one or more vertices). In other words, the spectrum of a 
Birkhoff integrable Hamiltonian system is obtained from the spectrum of a complete system 
by dropping part of the elements. 

Remark 2 The Dynkin type diagram determines only the angles between pairs of vectors 
in A. In order to reconstruct the ratios of lengths of vectors in A we assign to the ith 
vertex a coefficient proportional to the square of the length of Vj. This explains the numbers 
appearing on the vertices of the diagrams. 

We have to stress that this classification gives only necessary conditions for a system of 
type (1) to be Birkhoff integrable. The integrability for each system in the list should 
be established case by case. We give a brief history of the progress in this direction. As 
we already mentioned, the integrability of systems (a)-(g) was established in [1], [4]. The 
solution of these generalized periodic Toda lattices (associated with affine Lie algebras) was 
obtained by Goodman and Wallach in [14]. The graph (i) corresponds to a Hamiltonian 
system in two degrees of freedom with potential 

—*-) . 

The additional integral can be found in [19]. The integrability or non-integrability of systems 
(j) and (k) is still open. No Lax pair is known for either system. It is believed that system 
(j) is completely integrable, in fact integrability is established in [10] for the case n — 4. It 
is generally believed that system (k) in non-integrable. 

In this paper we deal with the integrability of system (h) which corresponds to the Hamil- 
tonian (8). Sklyanin in [27] indicated this system as another integrable generalization of the 
Toda lattice. The case n = 2 corresponds to the potential 

V = e qi ~ q2 + c x e 2q2 + c 2 e q2 + c^ qi + c 4 e~ 2qi . 

Annamalai and Tamizhmani [3] demonstrated the integrability of this particular case by 
using Noether's theorem. The second integral is of fourth degree in the momenta. 
The case n = 3 (as well as the general case) is treated in Ranada [26]. Ranada proved 
integrability by using a Lax pair approach. The additional integrals are of degree 4 and 6. 
In this paper we examine the integrability of the system 

n i n— 1 

H(q 1 ,...,q n ,p lj ...,p n )=Y / 7;P" + J2 e " +1 + ae ~ 2qi + ^ ■ ( 9 ) 

i=i z i=i 

Without loss of generality we may assume that a — f3 — 1. We note that if either a or 
(3 are zero then the system reduces to well-known cases of generalized Toda lattices.. The 
technique we use is entirely new and hopefully it will lead to further results for similar 
systems. The strategy is the following: We consider the D n Volterra system whose Lax 
pair formulation and multi-Hamiltonian structure was established recently in Kouzaris [18]. 
Using a procedure of Moser we transform the system into a new system in some intermediate 
variables (aj, bi): 

di = 2a 1 b 1 (10) 
hi = ai(bi-bi-i) i = 2,3, . . . ,n 
h n +i — — 2a n+ ib n 

k = 2(a 2 i+l -a]) i = l,2,...,n. 



We obtain a Lax pair for this system and the integrals of motion. The final step is the 
construction of a Flaschka transformation from the Hamiltonian system (9) to the system 
(10). The inverse of Flaschka's transformation provides for the necessary constants of motion 
in the variables (q,p) for the system (9). Thus, integrability is established. 
In section 2 we describe the construction of Bogoyavlensky-Volterra systems following [5], 
[6]. Bogoyavlensky constructed the systems using the root system of a simple Lie algebra and 
then through a change of variables (from q to Ui) he ended-up with homogeneous polynomial 
systems in the new variables U{. The construction of Lax pairs in the variables Ui is in [18]. 
In section 3 we present part of the results of [18], namely the D n case, since it is the only 
system needed for the purposes of the present paper. We have to point out that in [18] there 
is a complete treatment of the multi-Hamiltonian structure, Lax pairs and integrability of 
Bogoyavlensky-Volterra lattices. 

The results of section 4 are entirely new. They can be summarized as two transformations, 
one Moser-type from the Volterra D n lattice to system (10) and one Flaschka-type from the 
Sklyanin system (9) to system (10). Since Flaschka-type transformations are well-known, 
we describe briefly the Moser approach. 
Consider the system 

— ^ = Ui(u i+1 - Ui-i) i = l,...,n, (11) 
at 

where uq = u n+ \ = 0. This is the Volterra system, also known as the KM system and is 
related to the root system of a simple Lie algebra of type A n . The infinite KM-system was 
solved by Kac and Van Moerbeke [16] using a discrete version of inverse scattering. The Lax 
pair for system (11) can be found in [22]. The Lax matrix has the form 



ai 




V 



o 

a 2 
a 2 a 3 













\ 



a n _i J 



(12) 



where Ui = 2a 2 . Moser in [22] describes a relation between the KM system (11) and the 
non-periodic Toda lattice. The procedure is the following: Form L 2 which is not anymore 
a tridiagonal matrix but is similar to one. Let {ei, e 2 , . . . , e n } be the standard basis of R™, 
and E Q = {spane 2 j_i, i = 1,2, . . .}, E e = {spane 2i , i = 1,2, . . .}. Then L 2 leaves E Q , E e 
invariant and reduces to each of these spaces to a tri-diagonal symmetric Jacobi matrix. 
For example, if we omit all even columns and all even rows we obtain a tridiagonal Jacobi 
matrix and the entries of this new matrix define the transformation from the KM-system to 
the Toda lattice. We illustrate with a simple example where n = 5. 
In this case 



L = 



( o 

a i 



V o 



ai 

a 2 







a 2 

a 3 







a 3 

a4 



\ 





04 
0/ 



(13) 



and L 2 is the matrix 



( al 









a\ + a\ 


a 2 a 3 




aia 2 


a\ + a| 





a 2 a 3 


a§ + 04 






a 3 a 4 




a 4 



(14) 





a 3 a4 

Omitting even columns and even rows of L 2 we obtain the matrix 



aia 2 





a\a 2 
a\ + a\ 
a 3 a 4 





a 2 



(15) 



This is a tridiagonal Jacobi matrix. It is natural to define new variables A\ = aia 2 , A 2 



a 3 a 4 , Si 



< 1 • 



i? 2 = «2 + a 3> B 3 = aj. The new variables A 1 ,A 2 ,B 1 ,B 2 , B 3 satisfy the Toda 



lattice equations. 

This procedure shows that the KM-system and the Toda lattice are closely related: The 
explicit transformation which is due to Henon maps one system to the other. The mapping 
in the general case is given by 



Ai 



Bi 



\{u 2l 



- W21-2) 

The equations satisfied by the new variables Ai, Bi are given by: 



(16) 



B, 



Ai (B t , 

2 (A 2 - 



1 - Bi) 



These are precisely the Toda equations in [12]. We remark that the transformation (16) 
was first discovered by Henon. Henon never published the result but he communicated the 
formula in a letter to Flaschka in 1973. We refer to paper [7] for more details. In [7] one can 
find the multiple-Hamiltonian structure, higher Poisson structures and master symmetries 
for system (11). 

We would like to generalize the Henon correspondence (16) (using the recipe of Moser) from 
generalized Volterra to generalized Toda systems. The relation between the Volterra systems 
of type B n and C n and the corresponding Toda systems is in [8]. It is natural to attempt 
to find a similar correspondence between the Volterra lattice of type D n and the generalized 
Toda lattice of type D n . It is a surprising result, and this is the content of the present paper 
that the Volterra D n system corresponds not to the Toda D n system but to a special case 
of the Sklyanin lattice. 



2 Bogoyavlensky— Volterra systems 

Bogoyavlensky constructed integrable Hamiltonian systems connected with simple Lie alge- 
bras, generalizing the KM-system (11). For more details see ref. [5], [6]. In this section we 
summarize the construction of Bogoyavlensky. 

Let Q be a simple Lie algebra of rank n and n = {u>i j o; 2 , • • • , uj h } the Cartan-Weyl basis of 
the simple roots in Q. There exist unique, positive integers ki such that 



k uj + ki uji + 



• + k n u n = , 



where k = 1 and u is the minimal negative root. 
We consider the following Lax pairs: 



L = [B,L] , (17) 

n 

i=l l<i<j'<n 
i=l C « W 

Let Hbe a Cartan subalgebra of C? . For every root u; a G H* there is a unique H Wa e W such 
that co (h) — k (H Wa , h) V h E H, where is the Killing form. We also have an inner product 
on H* such that (u a , uj b ) = k (H Wa ,H u , b ). We set 



Cij 



1 if (oji,ujj) 7^ and i < j 
if (ui,cvj) = or i — j 
—1 if (uji,u>j) 7^ and i> j . 



The matrix equation (17) is equivalent to the dynamical system 



n he- 

c i = -Y t —- (18) 



We determine the skew-symmetric variables 



•Eij C-ijC^ Cj , Xij , 5 



which correspond to the edges of the Dynkin diagram for the Lie algebra Q, connecting the 
vertices uji and ujj. 

The dynamical system (18) in the variables takes the form 



x 



Xij ^ ' k s {xi s -\- Xjg) . (19) 



s=l 



We recall that the vertices u>i, ojj of the Dynkin diagram are joined by edges only if (u>i, uij) 7^ 
0. Hence x^ = if there are no edges connecting the vertices uji and ujj of the diagram. We 
call the equations (19) the Bogoyavlensky-Volterra system associated with Q (BV system 
for short). 

We shall now describe the BV system for each simple Lie algebra Q. The number of inde- 
pendent variables x^ (t) is equal to n — 1 and is one less than the number of variables Cj (t) . 
We use the standard numeration of vertices of the Dynkin diagram and define the variables 
life (t) = x^ (t) corresponding to the edges of the Dynkin diagram with increasing order of 
the vertices (i < j) . 

The phase space consists of variables Ui, with Ui > 0. In the following list we give explicit 
expressions for the Volterra systems in the variables Ui for each classical simple Lie algebra. 



A 



n+l 



O 



Ui 



UJ 2 

o 



U 2 



L0 3 

-o 



w = -(wi + wH W n+ i) 

fcj = 1 , i = 1, . . . , n + 1 






\i 




1 


j 


= i + l 


-1 


j 


= i-l 



Ui 



CiCi + l 



l,...,n 



Ui 



U\U 2 



Ui = Ui(u i+1 - Ui-i) 

Ki<n-\ 



B 



n+l 



O 



O 



«2 



W 3 



U r . 



o 



a; = - (a>i + 2u; 2 H h 2a; n +i) 

fci = 1, ki = 2 , i — 2,...,n + l 





1 j = i + 1 
-1 j = « - 1 



KI 



Ui Xi^i+i 



CiCi + l 



i = 1, . . . , n 



iii = ui (ui + 2u 2 ) 
u 2 = u 2 (2u 3 - ui) 
u n = 2u n ^\u n 
iti = 2ui(u i+1 - 
3 < i < n - 1 



a 



n+l 



O 



CJ 2 

o 



«2 



w 3 

-o 



o 



o 



uj = - (2o»i + 
&i = 2 , i — 1. 



-1 



1 



• + 2cU n + Wn+i) 

. , n , /c ra +i = 1 

J = 1 + 1 

j = i - 1 

= l,...,n 



ttx = 2u±u 2 

U n -1 = M n -l(«n - 2m„_; 

«« = -u n (u n + 2w„_ 

= 2Uj(lt i+ i - Ui-i) 

2 < z < n - 2 



n+l 



CJ 2 

o 



u 2 



u 3 

o 




ujq — - (ui + 2u; 2 H h 2o;„_i + w n + ^n+i) 

fci = 1 , £;„ = 1 , k n+ i = 1 , ki = 2,2<i<n — 1 
1 2 < j = z + 1 < n 
(i, j) = (n, n + 1) 

3 < i + 2 < j <n 

1 = (n- l,ra+l) 



CjCj + 1 



Cn-lC n +l 



ill = Mi (2w 2 + Mi) , U 2 = «2 (2^3 - «l) 

-Uj = 2-Uj(wj + i — 3 < z < n — 3 

Mn-2 = «n-2 + M„_l - 2« n _ 3 ) 

«n-l = «n-l (ti„ - 7J„_i - 2« n _ 2 ) 

«n = -«n («i> - «n-l + 2li„_ 2 ) . 



3 The Volterra D n system 

Consider the BV D n+ \ system (20) in the variables u k . We make a linear change of variables 

vi=u u v k = 2u k k = 2, . . . n - 2, v n ^ = w„_i, v n = u n , 
to obtain the equivalent system 



Vi = V! (vi + v 2 ) 

Vk = v k (v k+l - v k -i) k = 2, . 

V n -2 = V n - 2 (v n + V n -! ~ U„_ 3 ) 

i> n = -v n (v n - v n _ x + v n _ 2 ) . 



. n — 3 



(21) 



Before giving the Lax pair for the system (21) we introduce some matrix notations: 



Y k = 



~k 

iJvk~ 



£ / y/VkVk+l 

2 



O 












Y ° ~ 2 l -«i 



We also set 



-y/Vn-l iy/Vn-1 



Y=-l V V n-2Vn ^Vn-2Vn \ 

~ 2\ -y/V n -2V n ~l y/V n - 2 V n -l ) ' 

v n _i -v n \ 



W = - I " ra i 

2 \ ^ - v n -i o y 

Equations (21) can be written in a Lax pair form L = [L, B], where 
























X 












X 1 





X n -2 












X n -2 




















x 2 












x 2 






(22) 



-\\/V\V2 -\y/V!V 2 











Y 

















w 





Y n -3 











-Y* 






















— Yn-3 









Y 3 






















Y 2 












-Y 3 
























-Y 2 





Y 



We note that the entries of the first row and the first column of L and B are scalars while 
all the other entries are 2x2 matrices. 

The invariant polynomials of this system are given by the functions 

H 2 , ff 4 , . . . , H n -\ when n is odd, 

H 2 , H4, . . . , H n -i when n is even, 

where H k = \Tr(L k ) . 

We use the variables Cj, 1 < j < n + 1 of equations (18) in order to find a cubic bracket 7r 3 
of the BV D n+ i system. The dynamical system (18) in the case of the Lie algebra of type 
D n+ i can be written in Hamiltonian form dj = {cj, H}, with Hamiltonian 

n-l 

H = log ci + 2 lo S c j + lo S °n + log C n+ i , 

3=2 

and Poisson bracket 



{cj,c j+1 } = -{c j+1 ,Cj} = l j = 1,2,..., ra- 1 (23) 

{ c n-l) c n+l} — — { c n+li c n-l} = 1 

all other brackets are zero. In the new variables Vj (vi = c^c^ 1 , v k = 2cj~ 1 cj~+ 1 , 
k = 2, ... ,n — 2, v n _i = c^c" 1 , i> n = c^^c"^ ) the above skew-symmetric bracket, 
which we denote by 7r 3 , is given by 



{vi,v 2 } = v 1 v 2 (2v 1 +v 2 ) 

{vi,v i+ i} = ViV i+1 (vi + v i+1 ) , i = 2, . . . , n - 3 

{^n-2, l>n-l} = V n _ 2 V n _ x (2v n _ 1 + V n _ 2 ) 

{Vn-1, V n ] = 2v n -iV n (v n - 7J n _i) 

{^,^+2} = ViV i+1 v i+2 ,i = 1, ...,n- 3 

{i>„_2, V„} = W„-2Vn (Vn-2 + 2v n ) 

{v n _ 3 ,V n } = W„_3V„-2Vn ; 



(24) 



all other brackets are zero. As in the case of KM system we suppose that n is odd 
(n = 2m + 1) and we look for a bracket 7Ti which satisfies 7T 3 V H 2 = 7Ti Vi/4. 
We define 

n^2fc+l f . . 
for I < J , Tij = V 2i -1 , 

fc=i ^2fc 

and we define the bracket ni as follows: 

{vi,Vj} = (-l) <+j_1 T|^ +li [i+i] for l<i< j<n-2 , (25) 

(-l) i+n 

H,v n _i} = {^,d„} = r [i]+i,[2] for i = l,...,n-2 , 

{«n-l,Un} = -{Vn,«n-l} = ^(Un-V„_i) . 

To illustrate, we give the Poisson matrix of the bracket 7Ti in the case n = 7. 







r n 


-T12 


Tl2 


-Tl3 




5T13 


|ti 3 






-rn 


T22 


— T~22 


T23 




1 

~2 T 23 


1 

-2 T 23 




Tl2 — T22 





T22 


— T23 




1^ 
2 T 23 


1_ 

2 T 23 


7Ti = 




-Tl2 T 2 2 
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The following theorem is from [18]. 



Theorem 2 (i) 7Ti , n 3 are Poisson. 
(ii) The function 

n— 2 n— 2 

-H 2 = -Tr(L 4 ) = v n _ 2 v n + 2M n _if„ + + o XX 

4 » i=l 1 i=2 

is the Hamiltonian of the BV D n+ \ system with respect to the bracket 7Ti. 
(in) The function 

n-2 

F = (v n - u n _i) Y[ Vi , 

the Casimir of the BV D n+ i system in the bracket 

(iv) 7Ti , 7r 3 are compatible. 

(v) TT 3 VH 2 = TTiVi^ ■ 

4 From Volterra to Birkhoff 

We consider the Volterra D n+ \ system (21). We assume that n is odd, equal to 2m + 1 
and rename again the variables (i.e. use in place of v k ). We recall the equations for the 
system: 

iii = ui (ui + u 2 ) (26) 
ii k = u k (u k+1 - u k -i) 2 < k < n - 3 

U n -2 = «n-2 (U n + M n -1 ~ W„- 3 ) , 
«n-l = (u n - M„_i - W„_ 2 ) , 

U n = -U n (U n - M n _i + U n - 2 ) . 

We make the transformation (the analogue of Henon transformation (16) for the KM-system) 

ai = - (u n - u n -i) (27) 

aj = ^u n - 2j+2 u n „ 2j+1 j = 2,3,...,m 
i 

h = -^(u n + M„_l + u„_ 2 ) 

= -i(M n -2j+l + «n-2j) j — 2,3, ... ,171. 

This transformation is derived by mimicking the construction of Moser which takes the 
A n Volterra lattice to the Toda lattice, i.e. the construction that was described in the 
introduction. The formulas (27) are obtained by considering L 2 , where L is given by (22), 
and assigning suitable entries to the variables ai, bi. 
We calculate: 



(ll = - (U n - Un-l) 

= ^ l~ U n( u n ~ + U n - 2 ) ~ U n ^{u n - W n _i - U n _ 2 )] 

= ^(-«^ - «n-2«n + «Ll + M„-2« n -l) 

= 2a\b\. 

In a similar fashion we obtain the equations of motion in the new variables (a^, hi). 
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We have that L = [B, L] is equivalent to equations (28). 
The functions 

1 



H 2k = — tr L 



2 k 



2k 



k — 1, 2, . . . , m . 



are independent constants of motion. There is also a Casimir, which makes the system 
completely integrable, but these functions are enough for our purpose. 
We take H 2 as the Hamiltonian and define a Poisson bracket ix\ as follows: We consider 
the mapping from R 2m+1 to R 2m+1 given by (27). In other words it is the mapping which 



transforms the Ui to the new variables (aj,6j). The image of the bracket (25) (up to a 
constant multiple) is given by 



{aiA} = ai (30) 
RA} = -Oi i = 2,3,...,m 

{a i+1 ,bi} = ~a,i i = 1,2, . . .m - 1 

all other brackets are zero. The function 

C = a\a\a\ . . . a 2 rn a m+ \ 

is a Casimir. This is the Casimir F of theorem 2 in (a, b) coordinates. 

We also have involution of invariants, {Hi, Hj} = 0. We denote this bracket by n±. We have 

is equivalent to equations (28). 

Define a Hamiltonian system in R 2m with coordinates (qi, . . . , q m , Pi, ■ ■ ■ ,p m ) by 

1 m m—1 

H(q,p) = t EP? + E eqj ~ qj+1 + e ' 2qi + e29m - ( 31 ) 

2 3=1 3=1 

As in the case of the classical Toda lattice we make a Flaschka-type transformation 

ai = -^e-«\ a m+1 = -±=e^ (32) 
Oj = -e^*- 1- *) i = 2,3,...,m 

1 

&i = -gP* « = 1,2, ...,m . 

This is a mapping from R 2m — > R 2m+1 . Note that we are not using transformation (2) but 
a more traditional variation. 

We easily verify that Hamilton's equations for the variables (a,j,bj) are precisely equations 
(28). For example, 

1 OH 

ai = i^e~ qi qi = -cii — = -a x p x = -ai(-2&i) = 2ai&i . 

V2 <9pi 

We recall that the system (28) has H 2 , H4, . . . , H2 m as a set of integrals in involution. Re- 
verting back to the original variables pj) in R 2m we obtain m independent integrals in 
involution, and this proves the integrability of (31). For example 



#2 = E^ 2 + °? + 2 E a i +° 2 



m+l 



corresponds to the Hamiltonian (31). The integrals are of degrees 2,4,6, ... ,2m in the 
momenta. The Casimir C reduces to a constant equal to 22 m-i . 

We note that for the system (28) we may define a cubic bracket 7r 3 which satisfies the Lenard 
relation 

7T3 V H 2 = TTi V H A . (33) 

This bracket is the image of the bracket (24) under the mapping (27). The bracket 7r 3 is 
given by 

{(Xj, a i+ i} 
{cii, a i+ i} 
{bi,b i+1 } 
{ai,6i} 

RA} 

{ai, ^2} 
{di,b i+1 } 

{a 2 ,h} 
{a i+ i,bi} 

{^m+l, &m-l} 

all other brackets are zero. It is clear, by the way it was constructed, that this bracket is 
Poisson, compatible with 7Ti and that the integrals Hi are all in involution. We close with a 
few remarks: 

1. One could try to obtain a bi-Hamiltonian formulation of the system (28) following the 
recipe of [9]. The basic steps in the construction are the following: Write the second 
bracket 7r 3 in (q,p) coordinates (call it J 3 ) and define a recursion operator in (q,p) 
space by inverting the standard symplectic bracket (Ji). Define the negative recursion 
operator in (q,p) space by inverting the second bracket J 3 . Define a new rational 
bracket J_i by J_ x = JiJ^Ji- Finally, project the bracket J_i into the (a, b) space 
to obtain a rational bracket 7r_i. The result is a bi-Hamiltonian formulation of the 
system: 

7riV// 2 = 7T-lV# 4 • 

Note that a recursion operator in the (a, b) space does not exist since 7Ti and 7r 3 are 
both singular. So far, we were unable to compute the bracket J 3 and we are not certain 
that it can be computed. 

2. It is a straightforward recursive process to solve the Ui as functions of the at, bi using 
(27). After substitution of the values of a^, bi from (32) we obtain an expression of 
the Ui as functions of q iy pi. The formulas are too complicated but in principle the 
invariants, symmetries and higher Poisson brackets in the (q,p) space of (31) transfer 
to the corresponding ones for the Volterra system (27) via this mapping. 



= aia i+1 bi i = 2,3, . . . ,m - 1 

= 2aiCL i+ ibi i = 1 and m 

= 2a 2 i+1 (bi + b i+ i) i = l,2,...,m- 1 

= 2 ai {a\ + b\) 

= ai(a 2 i+b 2 i) i = 2, 3, . . . , m — 1 

= a m (a 2 m + b 2 m -a 2 m+1 ) 

= 2a\a\ 

= a 2 +1 cii i — 2, 3, . . . , m — 1 

= —a 2 (al + b\ — a\) 

= -a i+1 (a 2 +1 + b 2 ) i = 2,3, ...,m- 1 

= -2a m+ i(a^ +1 + 

= -a 2 +l a i+2 % = 1,2, . . . ,m - 2 



3. One may predict the degrees of the integrals by computing the Kovalevskaya exponents 
as in [11]. In our case the degrees of the invariants 2,4, .. . agree with the predicted 
values. 

4. The Casimir 

r i _ 2 2 2 

may be obtained in a different way using the following observation: The Casimir is 
the product of all the a« raised to certain exponents. We note that the exponents (not 
to be confused with the Kovalevskaya exponents) (1, 2, 2, . . . , 2, 1) can be determined 
from the condition 

Vi + 2v 2 + . . . + 2v7v_i + v^v = . 
This is always the case for systems of the form (1) if N > n. 

5. The procedure of this paper works equally well for the system 

H( qi , ...,q n , Pl , ... ,p n ) = J2 \p1 + E e ^ I+1 + ae ~ 291 + ■ (34) 

i=i 1 i=i 

The only difference is that we use the Volterra D n+1 system with n = 2m. 

6. The Hamiltonian (31) is positive, and the Casimir C keeps the variables from becom- 
ing zero. Thus, the energy surface is compact, and solutions lie on tori. This suggests 
that it should be possible to introduce a spectral parameter into the Lax pair and so 
get to Riemann surfaces and theta functions. 

Acknowledgments We thank the anonymous referee for useful remarks and corrections 
which made the paper more readable and precise. 
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